R Notebook

library(sf)
library(leaflet)
library(tidyverse)
library(rgdal)

Warning: package ‘rgdal’ was built under R version 4.1.1
Loading required package: sp
rgdal: version: 1.5-23, (SVN revision 1121)
Geospatial Data Abstraction Library extensions to R successfully loaded
Loaded GDAL runtime: GDAL 3.2.1, released 2020/12/29
Path to GDAL shared files: C:/Users/bewar/Documents/R/win-library/4.1/rgdal/gdal
GDAL binary built with GEOS: TRUE 
Loaded PROJ runtime: Rel. 7.2.1, January 1st, 2021, [PJ_VERSION: 721]
Path to PROJ shared files: C:/Users/bewar/Documents/R/win-library/4.1/rgdal/proj
PROJ CDN enabled: FALSE
Linking to sp version:1.4-5
To mute warnings of possible GDAL/OSR exportToProj4() degradation,
use options("rgdal_show_exportToProj4_warnings"="none") before loading rgdal.
Overwritten PROJ_LIB was C:/Users/bewar/Documents/R/win-library/4.1/rgdal/proj

st_transform(dz11, '+proj=longlat +datum=WGS84')

Simple feature collection with 6976 features and 9 fields
Geometry type: MULTIPOLYGON
Dimension:     XY
Bounding box:  xmin: -8.650007 ymin: 54.63324 xmax: -0.7246094 ymax: 60.86077
CRS:           +proj=longlat +datum=WGS84

lookup <- lookup %>% select(-2)

hb <- dz11 %>% 
  left_join(lookup, by = c("DataZone" = "DZ2011_Code"))

hb %>%
  filter(HB_Name == "Grampian") %>% 
  leaflet() %>% 
  addPolygons()

grampian <- hb %>% 
  group_by(HB_Name) %>%
  filter(HB_Name == "Grampian") %>%
  summarise(geometry = st_union(st_make_valid(geometry))) %>% 
  leaflet() %>%
  addPolygons() %>% 
  addTiles()

grampian

hb_polygons <- hb %>% 
  group_by(HB_Name) %>%
  summarise(geometry = st_union(st_make_valid(geometry))) %>% 
  ungroup()
  

hb_polygons %>% 
  leaflet() %>% 
  addPolygons()